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We study Einstein-Maxwell-dilaton theories with a cosmological constant and U{1)^ gauge symmetry, 
considering metrics asymptotically approaching the Lifshiftz metric. We study the dependence of the phase 
diagram on the value of the dynamical exponent. Along the way, we apply holographic renormalization and 
propose a counterterm valid for arbitrary dimension and dynamical exponent in our setup. 
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1 Context and setup 

This note is concerned with the study of black hole solutions to the Einstein-Maxwell-Dilaton system with 
negative cosmological constant, and in particular solutions in which the time and the space coordinates 
scale differently: t — > A^t, — > Ax* , where z is called the dynamical exponent. This kind of scaling is 
observed in many condensed matter systems, and has been investigated with great detail from a holographic 
perspective since the seminal work HI, in which one includes an extra holographic coordinate, scaling as 
r A~^r, with a boundary in which the field theory lives. The precise relation between the field theory of 
interest and the dual holographic theory is not yet known, and a construction of the holographic dictionary 
for this theories is under current investigation ||2l[3l|4|. 

In m, the authors proposed that gravitational theories realizing the Lifshitz symmetry should asymp- 
totically approach the so-called Lifshitz metric 

which has 1 + d{d + l)/2 Killing vectors, related to the generators of the symmetry algebra in the field 
theory (spatial rotations, spatial translations, time translations and dilatations) |l5l. Such a metric can be 
obtained in Einstein gravity only in the presence of matter fields. Most of the solutions constructed in 
the literature consider a system in which gravity is coupled to a massive vector field, described by the 
Proca action. If the mass of the vector field has the appropriate value the metric is asymptotically Lifshitz. 
The caveat with the Proca approach to asymptotically Lifshitz black holes is that only a few solutions 
are known analytically (and for specific values of dimensions, d, and the dynamical exponent z), and in 
top-down approaches the solution relies on numeric calculations ||6l . To find a generic solution for d > 2 
and z > 1 it has proven useful to consider instead the Einstein-Maxwell-Dilaton action with negative 
cosmological constant Q. 

Following the holographic philosophy, the construction of black hole setups asymptotically approaching 
([T]i are interesting to describe thermal field theories. Furthermore, if the effects of a global U (1) symmetry 
in the field theory are to be included, one has to consider a gauge field in the bulk of the d + 1-dimensional 
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gravitational theory. With this in mind we study the action IS) 



i?-2A- 1(90)2 



1 ^ 



Sgh 



(2) 



where the second term is the Gibbons-Hawking term, needed to have a well defined variational problem 
for the graviton. Below we will see that additional boundary terms must be added to obtain a finite on-shell 
action. These terms and the GH term do not affect the equations of motion derived from the action, which 
present a solution describing a charged black hole/brane 
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where we have set A = —{d + z — l){d + z — 2)/2^^ and (x is the amplitude of the scalar field. Notice 
that in the metric we have generalized the spatial section to admit d — 1-dimensional maximally symmetric 
manifolds with positive, null or negative curvatures (fc = 1, 0, —1 respectively). In the case fc = — 1 one 
of the gauge fields acquires a phase and we do not consider it further The case fc = can be obtained by 
taking the large-sphere-radius limit from the fc = 1 case ||9l- We will focus in the last case in the following, 
keeping fc explicit in some expressions. From now on we will set = 3 for concreteness. In the z ^ 1 
limit, in which ([T]) goes to the AdSd+i metric, this solution goes to the AdS-Reissner-Nordstrom solution. 
This model has been used to study the conditions that strongly interacting fermions should obey to satisfy 
the ARPES sum rules ifTol . 

The caveat this solution presents is the behavior of the Ai 3 gauge fields and the scalar <f>. These diverge 
at the boundary, and the variation relative to the Killing vector ~ ztdt — rdr + x^di is not vanishing, 
therefore spoiling the asymptotic Lifshitz symmetry algebra. We will not try to make sense of these fields 
in this note, and therefore the model we use must not be considered as a dual to a non-relativistic field 
theory, but as an IR effective theory that must be completed in the UV. 

From the solutions (|5]l and (|7]i we can define charges associated to the badly behaved radial electric 
fields (see equation ( fTSI ) below). These charges are fixed in terms of d, z and /i (contrary to the case in 
equation (|6]l, where it is given by a free parameter p2 = p). This suggests that these two fields should be 
considered in an ensemble in which their charges are kept fixed. The original Einstein-Maxwell-Dilaton 
action is appropriate for an ensemble in which the charges can vary ||9l, and therefore it must be modified. 
In other words, the boundary condition for the gauge fields imposed at the boundary of spacetime is not a 
Dirichlet boundary condition, but Neumann. To account for this we have to supplement the action with a 
boundary term to make the variational problem well defined. This is accounted by 
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where hmn is the boundary metric and n is a unit vector orthogonal to r = constant hypersurfaces. 

An alternative (equivalent) course of action is to perform a Legendre transform of the action with respect 
to the electric fields under consideration ifTTI . This transformation changes the equations of motion, but 
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(O-® is still a solution, only that the fields At,i and are not present in the new setup, since the 
solution has been already plugged into the Legendre transformed action. 



2 Holographic renormalization 

On-shell action As usual in holography, the evaluation of the on-shell action ^ diverges when one 
evaluates it at the boundary, and a renormalization procedure is needed. Since the publication of the 
solution ©-(ISll, several papers appeared dealing with the holographic renormalization of asymptotically 
Lifshitz spacetimes 12113] I?]. These, however, focus on the Proca action, and we will not consider them in 
the following. Indeed, the renormalization we describe in this section can be understood from the procedure 
for asymptotically AdS spacetimes with a few changes. 

To construct the needed counterterm action in asymptotically AdS spacetimes one has to solve a 
Hamilton-Jacobi equation ifTSl . which can be done systematically by expanding the new term in a deriva- 
tives series. For a counterterm depending only on gravitational terms this accounts to an expansion in 
powers of the boundary Ricci scalar contractions of the boundary Ricci tensor R{h) ,mnRJi^ , the 

boundary Riemann tensor, and generic combinations of them. Here we denote the boundary metric as 
hmn- The number of counterterms to add depends on the dimensionality of the theory under consideration. 
Generically, terms with i?"^^ (and combinations of a boundary Ricci scalars, Ricci tensors and Riemann 
tensors with all indices contracted) will be needed if the field theory dimension is d = 2 + 2a. 

In the present case such terms will regularize a field theory dimensionality d = 3 + 2a — z. This 
implies that to regularize the theory under consideration for generic z we need to consider an infinite 
number of counterterms. This is an impossible task, since for increasing value of a the combinations of 
curvature quantities increase in number (see flJl for the generic expression up to sixth order in derivatives 
in asymptotically AdS spacetimes). 

There is one property of our solution we can use to circumvent this, though. Focusing in the k — 1 case, 
the boundary metric has a Kt x S"^^^ topology and is static, implying that the boundary Riemann tensor 
takes a very simple form: any component with an index along the t direction vanishes, and the spatial 
part is that of a maximally symmetric d — 1-dimensional manifold with positive curvature. In particular, 
this means that any contraction of a boundary Riemann tensors, Ricci tensors and Ricci scalars will be 
proportional to R'^^i^y implying that the counterterm, expressed partially on-shell, is a power series in R^h)- 
For fixed d and z we can truncate the series and find what are the coefficients weighting each power 
of the boundary Ricci scalar to cancel divergences. When several examples are worked out, we can find a 
generic expression for the coefficient in front of every power of the boundary Ricci scalar, for arbitrary d 
and z. Assuming that these coefficients are the first coefficients of the infinite series in we induce the 
form of the 71*'* coefficient accompanying R?^^. in the series, and re-sum it. After this is done we find the 



This counterterm precisely cancels the value of the on-shell action (|9]l in the case where no black hole is 
present, meaning that we could have defined it by covariantizing the on-shell action (using he boundary 
Ricci scalar) in the zero temperature, zero charge setup. This is how the z = 1 case was found in ifTSl , 
and shows that the counterterm is valid in any number of dimensions and for arbitrary dynamical exponent 
z > 1, despite it being derived in the inductive way described before. 



This counterterm in the asymptotically AdS case, z = 1, was suggested in 1141 for d = 3, and in 1131 it was written for generic 
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For the thermal case we obtain a finite renormalized on-shell action, Sreii = Son-sheii + Sgh + Sc 
evaluated at a cutoff (and taken to the boundary) 



h 



One could have also considered the Legendre transformed action with respect to the physical charge p. 
It can be checked easily that the same counterterm does the job in this case, since the Legendre transform 
is equivalent to the inclusion of a boundary term (see equation ^ and the comments around it), which in 
our case give a finite contribution. After this consideration, the on-shell Legendre transformed action reads 

= TeCT^ ("■'^^" - - 4) - + . - - f/:f . (.2, 

Internal energy The trick used in the renormalization of the on-shell action is not useful when we try 
to renormalize the stress-energy tensor in the field theory. The reason is that we have used the fact that 
topology of the boundary metric is Mf x S"^^^ to express curvature invariants composed by the Riemann 
tensor, the Ricci tensor and the Ricci scalar as powers of the Ricci scalar alone. However, the calculation 
of the boundary stress energy tensor involves a variation of the action with respect to the boundary metric, 
and this result will be sensitive to the precise combination of curvature tensors that we have. Despite these 
considerations, let us write here the Brown- York tensor as derived by considering the counterterm ( fTOl i 



" 8^Gd+i 1^ ™" ^ (d-l)(d + z-3)2''™"j ■ ^^^^ 

In principle, T,„„ should have corrections involving terms proportional to contractions of curvature tensors 
with two free indices (for example R[h),mn or R(h),pmqnl^{h))- As an example, the lengthy expression one 
gets considering the most generic counterterm up to six derivatives of the boundary metric can be seen in 

m- 

Now, the staticity of our solution implies that any component of the Riemann tensor with an index in 
the time direction vanishes. This in turn means that the tt component can still be evaluated, since the 
corrections to ( fT3] l cancel for this component. With the evaluation of Tu we can calculate the renormalized 
energy density in the field theory straightforwardly, obtaining 

with dct hij involving only the spatial direction components of the boundary metric, ^ = 9t a Killing 
vector and k the unit vector orthogonal to t = constant surfaces. This result coincides withe the Komar 
mass result used in IB] by subtracting a reference background. 



3 Thermodynamics 

The temperature can be defined from the inverse of the euclidean time by the absence of singularities when 
performing a Wick rotation, giving an explicit expression T ~ T{rh, p, p) that can be found in |l8l|. The 
mass parameter 7ti > can be traded by the radius of the horizon, which is defined by hk{rh) = 0. The 
entropy is given, as usual, by the Bekenstein-Hawking formula S = Vd-irf^^^ / d+i- 

The charge carried by the gauge field A2 can be obtained from the Gauss law, and the chemical potential 
associated to this charge in the field theory is obtained, as usual in the holographic context, from the value 
of the corresponding gauge field at the boundary. These calculations give 
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where we have used the usual boundary condition A2,t(r/i) = 0. 

An analogous calculation for the chemical potential associated to the badly behaved fields Ai and A3 
implies a divergent value. However, as we work in an ensemble in which the charge associated to these 
fields is kept fixed, explicit factors of these chemical potentials will not appear in the thermodynamic 
relations. 

Thermodynamic potentials With the relations given above it is straightforward to check the relation 

W ^ TIren ^ E - T S - Q<^ , (16) 

with T the temperature. Actually, the differential thermodynamic relations in the grand-canonical ensemble 
(in which Q is allowed to vary and the potential $ is fixed) are also satisfied. We then conclude that the 
renormalized on-shell action ( fTTT ) (times the temperature) corresponds to the free energy of the field theory 
in the grand-canonical ensemble. 

For the Legendre transformed action (fT2l i one finds 

F = T{len - ire=n) ^AE^TS, (17) 

and looking at differential relations the identification with the free energy in the canonical ensemble is 
realized. This confirms the educated guess in IH, where these potentials were defined from the r.h.s. of the 
previous two equations, without going through the renormalization procedure. In ( fTTb we have subtracted 
the T = value to make contact with the analysis in that paper Here AE = E — Ef-xt is the energy with 
the value corresponding to the extremal solution subtracted. 

Phase diagrams We show the dependence of the phase diagrams with the dynamical exponent in figure 
[T] for both the canonical and grand-canonical ensembles. 

For 2 < 2 the situation is analogous to the AdS-Reissner-Nordstrom case discussed in ||9l (which 
corresponds to the z = 1 case in this note). There is a line of first order phase transitions in the grand- 
canonical ensemble separating the thermal Lifshitz solution (no black hole but periodic euclidean time) 
valid at low values of the temperature and the chemical potential, from the black hole setup. For large 
values of the chemical potential and vanishing temperature the dominating solution is a finite-entropy 
extremal black hole, which is expected to decay into an unspecified ground state. In the canonical ensemble 
there is a line of first order phase transitions between small and large black holes, terminating at a critical 
point above which we have a crossover. At Q ~ we have the analogous to the Hawking-Page transition 
lfT6l . The shadowed region in this diagram marks electrically unstable setups. 

For z = 2 the situation in the grand-canonical case stays unchanged, but in the canonical case the line 
of first order phase transitions gets reduced to just the critical point at Q = 0, i.e., to the Hawking-Page 
transition. The electric instability also disappears for this value of the dynamical exponent. 

Finally, for z > 2 the phase diagrams simplifies in both ensembles. In the grand-canonical one the 
first order transitions disappear and the black hole setup (extremal black hole setup in the T = case) 
dominates everywhere in the phase diagram, except at a critical value of the chemical potential (for T = 0) 
where the dominating solution is Lifshitz spacetime. In the canonical ensemble we do not find any phase 
transition. 

4 Conclusions 

We have studied the thermodynamics of the solution presented in IS) using holographic renormalization to 
calculate the free and internal energies of the system, related to the on-shell action and the tt component 
of the boundary Brown York tensor, respectively. 

The fact that a counterterm can be given for generic values of d and z is possible due to the extremely 
simple boundary topology of the solution, which allows to re-sum the infinite series of counterterms. The 
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most important limitation of this procedure being the impossibility to calculate the spatial components of 
the boundary stress-energy tensor 

In this renormalization procedure it is fundamental that we have performed the Legendre transformation, 
or equivalently the addition of the boundary term in canceling the divergent behavior of the gauge 
fields Ai and A3 in the on-shell action. We expect to revisit the role of the boundary-divergent fields and 
to provide an extended analysis of the renormalization procedure somewhere else. 
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